Surfatron acceleration along magnetic field by oblique electrostatic waves 
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Charged particles can be accelerated via surfatron mechanism along dc magnetic field by obliquely 
propagating electrostatic waves. In plasma, this mechanism can, in principle, produce an average 
parallel current, even when the wave frequency is much larger than the gyrofrequency and the wave 
phase velocity is much larger than particle initial velocities. 



I. INTRODUCTION 

It is well known that a charged particle can be ac- 
celerated resonantly by a wave along the phase front in 
the presence of weak dc magnetic field Bo [I]. This so- 
called surfatron acceleration was described in Ref. [2] for 
Bo perpendicular to the wave vector k. In particular, it 
was shown that speeds of the order of the E x Bo veloc- 
ity can be attained for nonrelativistic particles, E being 
the wave electric field. A more general case, with the 
angle between Bo and k being arbitrary, is more compli- 
cated, and it is only recently that the underlying physics 
was made transparent, namely, via Hamiltonian geomet- 
rical arguments in Ref. [3J. However, those arguments do 
not immediately apply to electrostatic waves, which also 
could drive surfatron acceleration in plasma (e.g., ion ac- 
celeration in the lower-hybrid frequency range). Hence, 
an operational understanding of the surfatron effect in 
oblique electrostatic waves is yet to be developed. 

Here, we study this effect in application to nonrela- 
tivistic particles. We proceed like in Ref. [3], but with- 
out using the Hamiltonian formalism explicitly, striving 
to render the qualitative picture even more transparent. 
For the particular case when k is perpendicular to Bo, 
the corresponding result of Ref. [2 is recovered. Oth- 
erwise, acceleration along magnetic field is possible, like 
in Ref. J3J; yet, the particle dynamics in an electrostatic 
wave can be qualitatively different from that in an elec- 
tromagnetic wave (Sec. III). We also find that, in plasma, 
surfatron mechanism can, in principle, produce an aver- 
age parallel current, even when the wave frequency is 
much larger than the gyrofrequency and the wave phase 
velocity is much larger than particle initial velocities. 

The presentation is organized as follows. In Sec. |Hj 
we introduce basic equations. In Sec. |III[ we identify 
and explain the possible types of particle trajectories. In 
Sec. IV we discuss current drive via surfatron accelera- 
tion. In SecfVl we summarize our main results. 



(For clarity, we do not constrain the wave to obey a spe- 
cific dispersion relation.) Then, the equation for the par- 
ticle momentum p reads as 



p = eV0 — fip x z 



(2) 



where f2 = eB /(mc) is the gyrofrequency, — e and m are 
the particle charge and mass, and c is the speed of light. 
Let us measure p in units mc, a in units mc 2 /e, Q in 
units u>, k in units oj/c, spatial coordinates in units c/u>, 



and time t in units 



Then, 



Px = -Pytt, 
Py = Px& + k y asimp, 
Pz = k z asimp, 

1p = 1 - kyPy - lt Z P Z , 

or, alternatively, 

p x + n 2 p x — —Qk y as'mip, (7) 
tjj + k 2 asinip = — ilk y p x . (8) 



(3) 
(4) 
(5) 
(6) 



Notice that Eqs. (|3j)-([8j) have an integral pi+p y + (p z — 
fcj 1 ) 2 — 2acosV> = const, owing to the energy conserva- 
tion in the frame where the wave is stationary [3] . Assum- 
ing that the wave is weak, the particle motion is thus ap- 
proximately bound to a sphere S p in p-space [Fig.[lja)], 



2,2,/ \2 2 

Px + Py + [Pz - U) = p , 



(9) 



with the center being at O s = (0, 0, u). Here p the radius 
determined by initial conditions, u — (cosa) -1 is the 
wave longitudinal phase velocity, and a is defined as 



sin a = k y /k, cos a = k z /k. 



(10) 



Below, we will assume, for clarity, that < a < tt/2, so, 
in particular, u > 0. 



II. BASIC EQUATIONS 

Suppose a dc magnetic field of the form B = z°B 
and a wave field E = — V<^, with some frequency w and 
the wavevector of the form k = y°k y + z°k z ; namely, 



III. SINGLE-PARTICLE MOTION 

For simplicity suppose k y ~ k z ; then, it is easy to see 
that the characteristic frequencies in Eqs. Q and (JsJ) are 
n and uje = k^fa. Of interest for us will be the essentially 
nonlinear regime, when 



4> = acosip, i(j — cut — k y y — k z z. 



(1) 



e = n/u E < i. 



(ii) 



FIG. 1: (Color online) Momentum space geometry, (a) Here S p is the sphere defined by Eq. (TsJ ; S is the resonance plane 
defined by Eq. ( |12[ ); C p is the circle obtained as the intersection of S p and S; L p is a Larmor orbit; Ri y 2 are the points where 
L p intersects with 3 (and C p also), (b) Here O s is the center of C p (and of S p also); the ring C is as defined in Sec. IV (c) The 
ring C with the radius p±o — p(p z o), m which case 3 is tangent to C p . 



[Notice that, at a <C p 2 and kp < 1, one also has « 1, 
so Eq. ( 11 ) automatically implies 51 <C 1 as well.] In this 



case, the particle nonlinear dynamics is quasiadiabatic 
and can be understood as follows. 



A. Two-scale oscillations 

First, let us consider the characteristic features of the 
two types of oscillations, employing Eqs. ^ and ijjj). For 
the motion with d t ~ w^, use that p x 3> £l 2 p x , so one 
obtains p x < flk y /k 2 ; hence, p 2 /a < 9 2 -C 1. The am- 
plitude of p y -oscillations is larger, namely, p y ~ p x /0 
[Eq. M]; yet still p^/a _ e" 2 (p|/a) < 1. Also, Eq. (j|) 
yields p z < ka/u>E for the quiver part of p z \ hence, 
p 2 /a < 1. Assuming the wave is weak, oscillations at fre- 
quencies of the order of loe are thus of small amplitude 
in p-space. In contrast, at dt ~ f2, we have V> wf^; 
then, at sin?/; ~ 1, one gets p x ~ p y ~ p z ~ V, where 
V = ka/VL is the characteristic electric-drift velocity [S]. 
Since V 2 /a = 6~ 2 2> 1, oscillations at frequencies of the 
order of can be of large amplitude in p-space. 

Since high-frequency oscillations of p x are negligible 
(albeit not so for p y and p z ), Eq. ^ yields that the 
dynamics in space 'F = (ip, tj)) is similar to that of a non- 
linear pendulum with a bias torque. Hence, two regimes 
are possible: (i) If p x > p c = k 2 a/(flk y ), there is no 
local equilibrium. Thus, all trajectories are untrapped; 
specifically, ip increases or decreases indefinitely, result- 
ing in that the right-hand side of Eq. ^ averages to 
zero on time scales of the particle gyromotion. Then the 
wave has little effect on the particle dynamics, which thus 
represents (weakly perturbed) Larmor rotation with fre- 
quency in the plane (p x ,p y ), with p z approximately 
fixed, (ii) If p x < Pc, most of trajectories are still un- 
trapped. However, if a particle is close enough to the 



resonance (i/> < oje), it can be trapped near one of the 
local equilibria, ip = 2im, where n is integer. [We assume 
a > 0; otherwise, the equilibria would be ip — 2i:{n + 1).] 
Since the parameters of the potential well for ip change on 
the time scale large compared to that of the bounce oscil- 
lations, the phase area J = <f ip dtp captured by the tra- 
jectory in ^ is an adiabatic invariant. Hence, the particle 
remains trapped as long as J stays less than the separa- 
trix value J a {p x {t)), or, equivalently, while p x (t) < p s (J), 
where p s (J) < p c is the function inverse to J s (j>x)- 



B. Types of trajectories 

Since ip of a trapped particle oscillates around a fixed 
value, averaging of Eq. (U over its fast motion yields that 
such a particle is bound to the "resonance plane" 2 in 
p-space [Fig. [TJa)], namely, 



1 - kyp y - k z p z = 0, 



(12) 



which passes through the center O s of each sphere S p . 
The intersection of S p and S is a circle C p with the same 
center O s and radius p, yielding that p x is bounded from 
above by p. Hence, if p < p s (J), the particle is trapped 
forever. If p s (J) < p < p c or p s {J) < p c < p, then de- 
trapping occurs at p x = p s {J). In particular, for most 
deeply trapped particles, the latter corresponds to de- 
trapping at p x — p s {J = 0) = p c . (This is the surfatron 
acceleration discussed in Ref. [2], where zero J and a 
are taken, also rendering p infinite.) Hence, the types of 
trajectories that can exist in the system are as follows. 

Type A. - Particles can be untrapped forever and 
rotate adiabatically in the plane (p x ,p y ) with frequency 
57; p z remains approximately fixed then [Fig. |2^a)]. For 
that, the particle Larmor orbit L p must not intersect S; 



3 




FIG. 2: Examples of particle trajectories of Type A, B, and C, respectively, obtained via numerical integration of Eqs. (|3|-(|6| 
with = 10~ 2 , k = 8, a = 70°, a = 5.3 x 1CP 5 , = 5.8 x 10 -4 , and u = 0.37. The initial momentum in Cartesian coordinates 
is (p_i_ o cos0, pj_osin0, p z0 ), where p z0 = 0, so p±(p z o) = ky 1 « 0.13. (a) p± = 0.9fcy 1 , 9 = 37r/2, V = 7r/8; (b) p± = fc" 1 , 
9 = tt/2, ^ = tt/8; ( c ) P-Lo = 2&" 1 , 9 = 3tt/8, ^ = 7vr/8. 



i.e., its radius p± must satisfy p± < p±(p z ), where 



P±(Pz 



1.-1 



1 - /c z p 2 



(13) 



Type 5. — Particles can be trapped forever and rotate 
adiabatically along circles C p in the plane S, specifically 
around the equilibrium at O s and in the direction deter- 
mined by the sign of fi [Fig. [2fb)] . The frequency of these 
oscillations in the limit of small J equals ficosa, as seen 
from Eqs. ffi and p| after linearization for ^ <C 1. 

Type C. — Particles can undergo nonadiabatic dynam- 
ics which is a mixture of trapped and untrapped motion 
[Fig. 2fc)]. Suppose a particle starts out on a Larmor cir- 
cle like of Type A yet with p± > p±(p z ). Then there are 
two intersections of L p and S, at some p x — ±p r , where 
the resonance condition ( 12 1 is satisfied; we will call those 



resonant points i?i,2- There, the particle crosses the sepa- 
ratrix in VP and thus has a probability of hopping on the 
trapped orbit C p , if arriving at the appropriate phase; 
namely, required is [cf. Eq. (18])] 



±p r w -pesin-^, 



(14) 



which also implies that p r < p c is needed. Then the par- 
ticle would have to continue along C p toward smaller \p x \ 
from, say, Ri and toward larger \p x \ from i?2- (Remem- 
ber that the direction of motion along C p is determined 
by the sign of only.) Yet notice that the phase area 
J that is captured at equals J s (p r ) [6-8J, and thus 
further increase of \p x \ would result in immediate un- 
trapping. Thus, actual trapping of an initially untrapped 
particle can occur at Ri but not at i?2- 

Hence, the overall corresponding motion consists of 
the gyromotion (at angular velocity Q) which is inter- 
rupted incidentally at R\; then the particle proceeds to 
i?2 along C p (at angular velocity of about f2 cos a) and 
gets untrapped again, returning to the same Larmor or- 
bit. Interestingly, for the type of waves considered in 



Ref. [3] realized is a very different scenario. Specifically, 
only large p x (rather than small p x , in our notation) are 
allowed there after trapping, thus resulting in particle 
jumping between two Larmor orbits with the same p±. 

Type D. — Particles can undergo nonadiabatic gyro- 
motion [p± > p±(p z )] with p r > p c . Then, trapping is 
impossible, but multiple scatterings off the resonance at 
p x = ±p r can cause diffusion along S p . Eventually, such 
trajectories can fill out the entire subset of S p that satis- 
fies p± > p±(p z ) [S]- This subset is further called S p and 
has a characteristic pitted-olive shape shown in Fig. [3] 



IV. AVERAGE CURRENT 

The asymmetry of S p along p z with respect to the ori- 
gin can be used to generate an average flow, or current, 




FIG. 3: (Color online) The pitted-olive-shaped surface S p ac- 
cessible for particles through diffusion along S p . 
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FIG. 4: Numerical calculation of (Ap z )(pxo) by tracing in- 
dividual trajectories of 2 16 particles (for each p±o), which 
were initially distributed in p-space over a ring C with fixed 
(pxo,Pzo) homogeneously in tp and 6. Shown is a snapshot 
at t = 2it/Q, for the same parameters as in Fig. [2] At 
p±o < PiPzo), C does not intersect E, so there is no current. 
At pj_o = p(Pzo), the current is maximal, because all trapped 
particles are of Type B. At p±o > p(Pzo), particles of Type C 
appear, reducing (Ap z ); yet, the latter still remains positive. 

parallel to B via surfatron acceleration of particles from 
an initially isotropic distribution. This is seen as follows. 

Ring distribution. - First, consider an ensemble of 
particles initially occupying a ring C in p-space, with 
some {p±o,PzO < u) and homogeneous distribution over 
ip and 9, where 9 is introduced via p X Q = p±o cos 9 and 
PyO = P±o sin 9. If p±o < p± {pzo) , then C does not inter- 
sect S, so all particles are of Type A and will approxi- 
mately retain their initial (p±,p z ), yielding no added cur- 
rent. Suppose now that p±o > p±(p z o), so £ intersects 
5 at some p x = ±p r , corresponding to # i:2 [Fig. [Tib)]. 
Then, trajectories of other types are realized, ana cur- 
rent can be generated. In particular, at p r > p Cl the 
current is produced through diffusion, for all trajectories 
are of Type D and thus eventually cover the entire S p , 
changing the average velocity by (Ap z ) > [5]. In con- 
trast, at p r < p c , the acceleration is ballistic and takes 
less time (of the order of Below, we assume this 

latter mechanism and discuss it in more detail. 

Clearly, for given p z0 , the maximum added current is 
generated when p±o = p±{p z o), i.e., when £ is tangent 
to 5, so 9\ — 9 2 [Fig. Iljc)]. This is because all trapped 
particles are either of Type B or remain on C forever, 
so Ap z > for each particle in this case. At larger 



p±o, some of trapped particles are of Type C. They will 
spend some time in the region p z < p z o, thus reducing 
the current somewhat. Still, even in this case one can 
show that (Ap z ) > 0, as also seen numerically (Fig. [4]). 

General case. — To assess the current that could be 
produced in plasma (i.e., by particles that initially oc- 
cupy a three-dimensional volume, rather than a ring, in 
p-space), one only needs to integrate over the contribu- 
tions yielded by each C(p±o,p z o)- Since each C produces 
(Ap z ) > 0, the total added current will also be nonzero. 
Remarkably, even particles with p z o <C u will contribute 
here (Fig. H). This makes surfatron acceleration very 
different from traditional current drive schemes [10] . be- 
cause, in the considered range (fi -C 1), those are bound 
to rely on Cherenkov resonance. 

Finally, notice that electromagnetic waves of the type 
discussed in Ref. [3] can also drive current similarly. In 
fact, since no particles attain negative Ap z in that case 
(assuming p z0 < u and except in the diffusive regime), 
the total added current may be even larger. However, the 
actual amount of current that can be produced through 
surfatron acceleration in plasma is yet to be calculated. 



V. CONCLUSIONS 

In this brief note, we study surfatron acceleration of 
nonrelativistic particles by electrostatic waves propagat- 
ing obliquely to weak dc magnetic field B . For the par- 
ticular case when the wave vector k is perpendicular to 
Bo, the corresponding result of Ref. [2] is recovered. Oth- 
erwise, acceleration along magnetic field is possible, like 
in Ref. [3]; yet, the particle dynamics in an electrostatic 
wave can be qualitatively different from that in an elec- 
tromagnetic wave. We also find that, in plasma, surfa- 
tron mechanism can produce an average parallel flow, or 
current, even when the wave frequency is much larger 
than the gyrofrequency and the wave phase velocity is 
much larger than particle initial velocities. 
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